Evidence for a Peierls phase-transition in a three-dimensional multiple charge-density 

waves solid. 
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The effect of dimensionality on materials properties has become strikingly evident with the recent 
discovery of graphene. Charge ordering phenomena can be induced in one dimension by periodic 
distortions of a material's crystal structure, termed Peierls ordering transition. Charge-density waves 
can also be induced in solids by strong Coulomb repulsion between carriers, and at the extreme limit, 
Wigner predicted that crystallization itself can be induced in an electrons gas in free space close to 
the absolute zero of temperature. Similar phenomena are observed also in higher dimensions, but 
the microscopic description of the corresponding phase transition is often controversial, and remains 
an open field of research for fundamental physics. Here, we photoinduce the melting of the charge 
ordering in a complex three-dimensional solid and monitor the consequent charge redistribution by 
probing the optical response over a broad spectral range with ultrashort laser pulses. Although 
the photoinduced electronic temperature far exceeds the critical value, the charge-density wave is 
preserved until the lattice is sufficiently distorted to induce the phase transition. Combining this 
result with ab initio electronic structure calculations, we identified the Peierls origin of multiple 
charge-density waves in a three-dimensional system for the first time. 
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Charge ordering phenomena occurring upon symme- 
try breaking are important in solids as they give rise to 
new current and spin flow patterns in promising materials 
such as organic conductors [I], multilayered graphene [2] 
and transition metal oxides j3]. The possibility to in- 
vestigate the microscopic steps through which such or- 
dering transition occurs gives also the opportunity to 
speculate on more general aspects of critical phenom- 
ena. Charge-density waves (CDWs) [H03, sandpile au- 
tomata [6] and Josephson arrays [7] have been investi- 
gated in relation to the scale invariance of self-organized 
critical phenomena |8{ , of which avalanches are dramatic 
manifestations [9 . In one dimension, Peierls demon- 
strated that at low temperature an instability can be 
induced by the coupling between carriers and a periodic 
lattice distortion. Such an instability triggers a charge 
ordering phenomenon and a metal-insulator phase tran- 
sition, called Peierls transition, occurs [10 . Like for 
Bardeen-Cooper-Schrieffer (BCS) superconductors, such 
an electron-phonon interaction driven transition is ex- 
pected to be second order [10]. While this situation is 
fairly established in mono-dimensional organic materials 



|11) , increased hybridization leading to higher dimension- 
ality of a solid perturbs this scenario and makes the as- 
sessment of the microscopic origin of charge localization 
phenomena more difficult pJJ-flTj. 

Contrary to other low-dimensional CDW [TS] systems 
studied so far by time-resolved spectroscopies [TMT9"] . 
Lu5lr4Siio presents a complex three-dimensional struc- 
ture with several substructures such as one-dimensional 
Lu chains and three-dimensional cages in which a vari- 
ety of many-body effects (including superconductivity) 
originate j5U] (Fig. [6]A). Although a CDW occurs below 
Tcdw= 83 K, its microscopic origin is still debated since 
this transition is first order and isotropic without hystere- 
sis [51 [T7] , in contrast to the standard Peierls paradigm. 
Also, owing to their complexity, a detailed microscopic 
description of the properties of these solids is still lack- 
ing and their band structure has never been reported. 

The charge redistribution induced by a phase transi- 
tion in a solid can be obtained via the optical frequency- 
sum (f-sum) rule [6j, which states that the integral over 
all frequencies of the optical conductivity (cti(cl>)), termed 
spectral weight (SW), is constant and corresponds to the 
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number of electrons per unit cell: 

SW = ^N eff = j~v^) d « (1) 

where m e is the free electron mass, V the unit cell vol- 
ume and N e ff the total number of carriers. Performing 
the integral up to a frequency cut-off tt c (partial sum- 
rule) gives the number of carriers contained in the elec- 
tronic levels included in this energy range [22]. Opti- 
cal weight redistributions are known to accompany every 
BCS-like phase transition, as described by the Ferrel- 
Glover-Tinkham sum-rule [53]. Recently, it has been 
shown that information on the SW can be obtained in a 
model-independent fashion based on the analytical con- 
tinuation of holomorphic functions [4] . 

In this study we aim at obtaining the temporal evo- 
lution of the partial sum rule during the photoinduced 
annihilation of the CDW order parameter to identify the 
different states through which the electronic and the lat- 
tice structure evolve. To achieve this, the changes in 



the optical constants are monitored over a broad en- 
ergy range with a temporal resolution better than the 
electron-phonon coupling time. We performed pump- 
probe reflectivity using a broad (1.5-3 eV) supercontin- 
uum of 50 fs polarized pulses as a probe [33] ■ The tem- 
perature of the sample was controlled between 10 K and 
room temperature. The 1.55 eV pump fluence was var- 
ied between 0.8 and 3.1 mJ/cm 2 , corresponding to an 
absorbed fluence between 0.4 and 1.55 mJ/cm 2 ; in these 
conditions and for an initial temperature of 10 K, the 
electronic temperature reaches values between 302 and 
1100 K (see the transient temperature analysis below), 
corresponding to energies in the range 26-94 meV. There- 
fore, for the lowest pump fluences, the absorbed energy 
lies below the CDW gap estimated by static optical spec- 
troscopy (80 meV) [BJ. Our ultrafast reflectivity data 
were combined with static optical spectroscopy [BJ and 
new ah initio band structure calculations, allowing the 
identification of the spectroscopic signature of the CDW 
melting, the description of the orbitals and atomic mo- 
tions involved and their temporal evolution. 
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FIG. 1: (A) Experimental geometry and crystal structure; the pump and probe polarizations can be varied independently. (B) 
Static reflectivity measurements; vertical bars indicate the optical absorptions of the solid, the vertical dashed line the energy 
of our photoexciting pulse, and the probe supercontinuum range is highlighted as a light-blue shaded area, covering a spectral 
region of large anisotropy. (C) Calculated reflectivity along the a-axis and the c-axis. {D-G) Time-resolved reflectivity images 
for a pump fluence of 3.1 mJ/cm 2 and (D-E) T = 10 K, (F-G) T = 290 K. White traces in (D-G) are profiles extracted at 700 
fs, whose baselines are shown as dashed red lines. 
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The static optical reflectivity of Lu5lr 4 Siio for light polarized along the one-dimensional chains (c-axis), and 
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perpendicular to them (a-axis) is displayed in Fig. |6]S. 
Fig. |6]C shows the calculated reflectivity which agrees 
with the experiments; some discrepancies are observed 
at low energy along the c-axis where correlation effects 
are expected due to the presence of the monodimcnsional 
chains. Figs [6](X>-GJ show the transient optical reflec- 
tivity for different temperatures and light polarizations. 
The low temperature dynamic reflectivity decays faster 



than at room temperature (200 fs compared with 300 fs, 
for the same pump fluence and energy of 1.9 eV). Its spec- 
tral behavior is also rather different: at 10 K, a large re- 
gion is observed in which the photoinduced change in re- 
flectivity switches sign (see white traces) . The timescales 
and the absolute magnitude of the changes also depend 
on the polarization of the probe. 
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FIG. 2: Spectral Weight analysis. (A-C) 3D color maps of time-resolved N e ff. (D-E) Difference in N e ff between the normal 
(T — 120 K) and the Charge-density wave (T = 60 K) states for the two orientations, and between two different temperatures 
in the normal state (290 K - 120 K) (solid lines); time-resolved iV e // for T = 10 and 290 K, and probe parallel to c or a-axis 
(circles), at time delays of fs (D) and 1 ps (E). 



The evolution of N e ff is obtained from the transient 
reflectivity both via a Drude-Lorentz model (3D color 
maps in Fig. [8]) and the model-independent approach 
(white lines with error bars in Fig. |8]A and B) [53]. In 
these plots, the integral in Eq. [T] is displayed as a func- 
tion of the frequency cut-off, emphasizing the spectral 
regions of positive and negative change of the effective 
number of carriers. When the sample temperature is 
below Tqdw an d the probe light is polarized along the 



c-axis (Fig. MB), the N e ff change, defined as (iV e //(t,w)- 
N e ff(t< 0, uj)), becomes negative after 200 fs. Instead, it 
is always positive when the polarization is perpendicular 
to the c-axis (Fig.[8]A), or at room temperature (Fig.[8]C). 
This sign change at low temperature along the c-axis sig- 
nals a SW transfer from energies higher than 1.5 eV, 
typical of strongly correlated solids [26] . In Fig. [8]D and 
E, the quantities {N ef f (T 2 ,uj)-N eJ f {T u u))/N ef f(Ti,w) 
are displayed for T x = 60 K and T 2 = 120 K, with light 
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polarized along the c-axis (red line) or the a-axis (black 
line), and T x = 120 K and T 2 = 290 K for light polarized 
along the c-axis (light blue line) (data taken from [6]). 
The SW transfer across the CDW temperature along the 
c-axis becomes negative above 1.5 eV, as observed in the 
transient data; therefore, the spectroscopic signature of 
the CDW melting is the redistribution of carriers above 



1.5 eV where bands composed of Lu and Ir 5d states are 
found. Remarkably, we observe that this feature appears 
several tens of fs after laser excitation, suggesting that be- 
fore the CDW melts, an unconventional state of matter 
is produced in which some charge order survives despite 
the large energy deposited in the electronic subsystem. 




FIG. 3: (A) Time-dependent lattice (black) and electronic (blue) temperatures, obtained by solving the three-temperature 
model for T = 10 K, F = 3.1 mj/cm 2 and an energy of 1.84 eV, and corresponding transient Neff. The dashed line represents 
the zero baseline for N e fj and the CDW transition temperature. (B-C) Fluence dependence of the transient reflectivity traces 
at 10 K and 1.77 eV (B) and 2.21 eV (C) probing wavelengths. (D) Temperature dependence of the transient reflectivity at 
3.1 mj/cm 2 pumping fluence and 2.21 eV probing wavelength. 



To verify that indeed the behavior of the optical con- 
stants in the 2-3 eV region is mainly sensitive to the CDW 
melting, we measured the transient reflectivity spectra at 
different fluences (0.8 to 3.1 mj/cm 2 ) and temperatures 
(10 K, 80 K and 300 K). In Fig. [7]B and C, the transient 
reflectivity at 1.77 and 2.21 eV are shown for the different 
pump fluences. At 1.77 eV, the reflectivity changes are 
always positive, while at 2.21 eV a negative dip develops 
at fluences higher than 2.1 mj/cm 2 . In Fig.[7]D, the tem- 
perature dependence of the transient reflectivity is also 
shown, at a probe energy of 2.21 eV and for a fluence of 
3.1 mJ/cm 2 . Also in this case, the change of reflectivity 
has a negative part, associated to the CDW melting, at 
10 K, while it is always positive at temperatures higher 
than T C dw- 

The time scales involved in the melting of the CDW 
are obtained through a transient temperature analysis 



using the three-temperature model H3], appropriate 
for anisotropic materials with selective electron-phonon 
coupling [TTJ [TS1 OUJ . This model describes the energy 
transfer from excited carriers to the lattice as an effective 
heat exchange between them. We obtained the tempo- 
ral evolution of electronic and lattice temperatures |23j 
and display them in Fig. [7]/4 together with the transient 
AN e ff at 10 K and polarization along the c-axis. 

The CDW melting, indicated by AiV e / f becoming neg- 
ative at high energy, occurs after the lattice reaches a 
temperature above Tqdw- While the electronic tem- 
perature T e crosses Tcdw and reaches values as high 
as 1100 K, the system exists in a transient state where 
the signature of the molten CDW state is still absent. 
These simulations provide an estimate of the Eliashberg 
electron-phonon coupling parameter A ~ 0.6 and a frac- 
tion of strongly coupled modes ~ 20%, in agreement with 
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the value of A = 0.66 determined by combined heat ca- 
pacity and magnetic susceptibility measurements under 
pressure [10]. The rise time of our AR/R signal is in our 
temporal resolution range, i.e. 40-50 fs. This, together 
with the three-temperature model simulations, provides 
an estimate of the excited electron thermalization time 



around 50 fs, much faster than the CDW melting time. 
This rules out any possible electronic melting as in this 
case a clear signature of the molten CDW state should 
be present as early as the electron thermalization time, 
like recently observed in TaS2 [55] . 




FIG. 4: Transient state. (A) Transient reflectivity spectra at 60, 120, 200 and 500 fs respectively, colored asterisks. Static 
temperature difference spectra, R(150 K) - R(25 K), continuous line. (B-D) Cartoon of the CDW melting dynamics. The 
colored areas indicate the regions where excitation perturbs the charge order. 



To speculate on the microscopic nature of the interme- 
diate state, we consider that during the first 200 fs the 
electronic structure is relaxing after being put abruptly 
out-of-equilibrium by the pump pulse. By comparing 
the spectra of the photoinduced reflectivity changes at 
different time delays to the difference spectra between 
the static reflectivity recorded at different temperatures, 
the impact of the out-of-equilibrium electronic structure 
can be estimated. In Fig. 4^4, the transient reflectivity 
spectra at 60, 120, 200 and 500 fs after excitation are 
displayed together with R(150 K) - R(25 K). When the 
electronic structure and the lattice are at equilibrium, the 
temperature difference spectrum exhibits a positive peak 
around 2 eV and a minimum around 2.7 eV across the 
phase transition (red line). This behavior is quite similar 



to that of the transient reflectivity after 500 fs, at which 
time the electronic structure is likely at equilibrium with 
the lattice. AR/R shows also a positive peak around 2 eV 
and a negative minimun around 2.7 eV; this similarity in- 
dicates that photoexcitation and temperature induce the 
same carriers redistribution in the solid. At early times, 
while the shape of the transient spectra is still similar to 
the static one showing a positive peak around 2 eV and 
a minimum around 2.7 eV, a quantitative discrepancy is 
more pronounced and could originate from the very large 
electronic temperature jump induced by laser excitation 
(> 1000 K) compared to the static case in which the sys- 
tem is heated by a hundred degrees. These spectra sug- 
gest that, consistently with the three-temperature model 
analysis, a thermal electronic distribution of carriers is 
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established within few tens of fs. 

It is important to note at this point that while a neg- 
ative reflectivity change can be associated to the phase 
transition, the complete melting of the CDW state can 
only be claimed when the optical SW, which is the cumu- 
lative effect of the positive and negative changes, becomes 
negative. At a fluence of 3.1 mJ/cm 2 , around 1 carrier 
per unit cell is photoexcited. Photoexcited carriers create 



"hot areas" in which the CDW can be perturbed. This 
electronic perturbation diffuses in space via its interac- 
tion with the lattice and after some time, a thermally 
CDW molten state is reached. The delay of this process 
is governed by the microscopic details of the coupling 
between carriers and phonons in the CDW state. This 
scenario is pictorially represented in the cartoon in Fig. 
AB-D. 
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FIG. 5: Electronic structure calculations. (A) 3D map of the Fermi surface. (B) Band structure with atomic character, and 
partial density of states. (C) Representation of the main atomic orbitals (Ir 5d z 2 and Lu 5d yz ) forming one of the nested bands, 
and atomic movement of two of the most coupled phonon modes (thick orange and yellow arrows). Lul atoms form chains 
along the c-axis, Lu2 the ring in the (a, b) plane and Lu3 (not represented) are more three-dimensionally distributed. 



The presented situation is very distinct from an elec- that case the negative optical SW should be observable 
tronically driven melting of the charge order, since in without any delay; therefore, even though an instanta- 
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neous (< 60 fs) CDW melting would locally happen upon 
light excitation, it should propagate and melt the whole 
charge order within the electronic thermalization time, 
which we found to be around 60 fs. Such a situation is 
observed for example in TaS2 |33| . where a non-Peierls 
phase transition is claimed to take place. 

Further microscopic details are obtained through ab 
initio electronic structure calculations [23] , performed via 
both an all-electron full-potential linearized augmented 
plane- wave method [IS] and via pseudopotentials [IB] , 
Eleven bands having mixed character (color coded in 
Fig. [9]B) cross the Fermi level. The density of states 
(DOS) (Fig.[9]B) shows that the Fermi surface has mostly 
Ir 5d character. Some Lu character is found, coming from 
the bd orbitals strongly hybridized with both Si and Ir, 
producing delocalized electrons as observed in transition 
metal-silicides 36 . The Fermi surface is presented in 
Fig. [9]A, and from calculations of the generalized suscep- 
tibility, we identified two bands having nesting vectors of 
2/7 and 3/7 of c*, as observed by X-ray diffraction below 
Tcdw 53- This suggests the presence of multiple CDWs 
occurring at the same transition temperature. These 
two nested bands account for approximately 30% of the 
Fermi energy DOS, in agreement with estimates based 
on magnetic susceptibility measurements jlO] . While one 
of them comes mostly from Ir bd z 2 and Lul 5d yz (Lul 
atoms being along the chains [T7] , Fig. [9]C) , the second 
one is mainly formed by Lu2 5d yz -5d z 2 and all the 5d 
orbitals of both Lu3 and Ir. This complex and three- 
dimensional nature of the orbitals reveals the possibility 
of a 3D Peierls transition. We also show in Fig. [9]C the 
motion of two stronglycoupled optical phonons, yellow 
and orange arrows respectively, having frequencies (1.4 
THz and 2.8 THz) close to the coherent oscillations ob- 
served in the same material [37]. Ab initio calculations 
also estimate the electron-phonon coupling parameter of 
each individual lattice mode, yielding a total A s» 0.6 for 
the 18 % of the most coupled modes, in agreement with 
the three-temperature model simulations. 

The orbital occupation-number change induced by the 
CDW melting is signalled by the optical SW transfer 
around a pivot energy close to 1 eV, both in static 
and time-resolved data. Upon excitation, acting as a 
photodoping process, carriers are transferred from the Lu 
and Ir 5d states forming the CDW [23] . This transfer de- 
populates the nested bands and redistributes the excited 
carriers above the Fermi level, increasing the low-energy 
SW as well as the effective electronic temperature. The 
observation of a delay between photoexcitation and the 
melting of the CDW strongly suggests the Peierls nature 
of the phase transition in these compounds, rather than 
a pure electronic effects as in different systems [38]. In 
fact, while CDW are known to couple to specific lattice 
distortions, the intimate microscopic origin of the mech- 
anism described in [10| is that a periodic modulation of 
the lattice induces the Fermi surface instability causing 



the consequent charge ordering. In other words, to un- 
derstand correctly the origin of such a phase transition, 
a notion of causality must be invoked, which is directly 
obtainable via time-resolved experiments. Indeed, one 
must wait until the atoms have moved sufficiently to sup- 
press the distortion promoting the CDW state [TH [3J5] ; 
the latency time can be of the order of one quarter of the 
corresponding phonon period. In our case, 200 fs gives 
a frequency of 1.4 THz, corresponding to the mode rep- 
resented in orange arrows in Fig. [9]C. While the direct 
involvement of the lattice in the CDW formation reveals 
the Peierls nature of the metal-insulator transition in this 
material, its weakly one-dimensional electronic structure, 
together with the signatures of a first order CDW tran- 
sition and the suggested presence of multiple CDWs, are 
at odds with the conventional scenario describing Peierls 
transitions |15j . Recently, the notion that Peierls tran- 
sitions can happen in higher dimensions has been chal- 
lenged by new ultrafast experiments in TaS2 [33 , which 
has a 2D charge ordering transition. Also in this case, 
the delay between excitation and melting has been used 
as the key-argument for determining the microscopic na- 
ture of the phase transition. Our experiments contribute 
to this exciting and current debate by showing evidence 
for a Peierls transition in a three-dimensional system. 

The proper description of phase transitions is impor- 
tant for fundamental physics; in fact, phase transitions 
in complex systems often evolve through several uncon- 
ventional intermediate states |40j . Critical phenomena 
involving abrupt perturbation of ordered systems are ob- 
served in giant land-slides [5] and sandpile automata [BJ. 
In these cases and in CDWs, an external perturbation 
may induce a transient state in which the order is pre- 
served while undergoing a collective reorganization [5] . 
The speed of rock-slides, CDW slides, and flux jumps 
in the vortex state of type II superconductors, are self- 
organized critical phenomena and obey similar critical 
laws [U [7] [HI E]. In our experiments, the disappearance 
of the charge order manifests itself after a significant de- 
lay from the abrupt perturbation, and in the intermediate 
state the order seems preserved despite the large energy 
deposited. 
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EXPERIMENTAL SETUP 

We performed an experimental study of the transient 
optical properties of Lush^Siio using pump-probe broad- 
band reflectivity. Our experimental setup, located at the 
Laboratory of Ultrafast Spectroscopy, EPFL, Switzer- 
land, is presented in Fig. SI. An amplified Tksapphire 
laser system provides 1.55 eV pulses of 40 fs duration at a 
repetition rate of 1kHz. This laser beam is separated into 
two paths by a beamsplitter. The pump beam, whose 
energy is fixed at 1.55 eV, and the probe beam which 
is converted into a supercontinuum by a CaF2 nonlinear 
crystal. The probe spectrum covers a range between 1.5 
and 3.2 eV. The fluence of the white light probe is several 
orders of magnitude lower than the pump beam and has 
no effect on the measurements. The pump pulse is used 
to excite the system at the time t = 0, and the conse- 
quent variations of reflectivity are recorded with a delay 
t by changing the optical path difference between pump 
and probe (which can be achieved by setting a variable 
delay in the probe arm). 

Small reflectivity changes (down to the order of 10~ 3 ) 
are detectable owing to the signal-to- noise ratio of 10 4 , 
obtained by the use of an optical chopper on the pump 
beam which reduces its repetition rate by a factor of two. 
This allows the recording, for each time delay, of a single- 
shot spectrum with and without the pumping pulse. Fur- 
thermore, the variations of the white light intensity are 
taken into account using a reference signal, measured just 
before the sample position. 

Two identical spectrometers are used to detect the 
transient reflectivity signal and the reference one. By 
dispersing the incident light with optical gratings and 
sending it to an array detector, we obtain the complete 
spectrum for each single pulse with an energy resolution 
better than 1 meV. Subsequently, the group velocity dis- 
persion introduced by the chromatic aberration of the 
optics is corrected before the transient reflectivity anal- 
ysis. 

The pump-beam diameter on the sample surface was 
250 /im, while that of the probe was set to 100 /im in 
order to probe a uniformly excited area. Both spot sizes 
and spatial overlap were rigorously verified using a beam 
profiler and a charge-coupled device camera. During this 
experiment, the pump fluence was varied between 0.8 and 
3.1 mJ/cm 2 . 

Prior to these experiments, the sample was oriented us- 
ing X-ray diffraction and polished to obtain good optical 
quality. Its surface contained both (001) and (100) di- 
rections. Therefore, pump and probe polarizations were 
changed independently with quarter- wavcplates and were 
oriented either parallel or orthogonal to the c-axis (con- 
taining the Lu chains). Finally, in order to cool down to 
10 K and precisely monitor the sample temperature, the 
latter was placed in a liquid He cryostat. 



Such kind of experimental setup has been successfully 
employed to investigate photoexcited states in strongly 
correlated systems, such as the charge-transfer insulator 
CuGeOs [SI], as well as in bulk graphite and graphene 
paper [S2]. 

SPECTRAL WEIGHT ANALYSIS 

The definition of the optical spectral weight (SW) in 
a solid, given in Eq. (1) of the main article, makes 
clear that the evaluation of N e ff requires the knowl- 
edge of the real part of the optical conductivity down 
to zero frequency. Since this is not experimentally fea- 
sible, common procedures (including extrapolations or 
modeling based on the Drude-Lorentz formalism) are in- 
voked to obtain the whole <j\[oj) spectrum. Recently, 
it has been shown mathematically that according to the 
theorems on analytical continuation of holomorphic func- 
tions, and taking advantage of the Kramers-Kronig re- 
lation between the real and the imaginary part of the 
complex optical conductivity (c(uj) = <j\(oj) + ia-iioS)), 
the partial SW integral can be re-written as a function 
of a limited spectral interval in the form of Carleman- 
Goluzin-Krylov equations [S3]: 

rui max 

SW(n c ) = Re { lim / uju', fl c )a(uj')duj'} (2) 

where ii n (o/,Q c ) = J Q C Q n (tj', u>)du> are properly cho- 
sen kernel functions. In the case of real experimental 
data (with finite energy resolution and statistical noise), 
Eq. ([2]) cannot be strictly applied, and a numerical pro- 
cedure based on the idea of the analytical continuation 
has to be invoked in order to obtain the SW integral in a 
model- independent fashion [S4]. A code called "Devin" 
has been developed which allows the SW integral to be 
extracted from the complex dielectric function in a lim- 
ited data range in a model-independent fashion [S4] . 

We performed a Drude-Lorentz analysis of our exper- 
imental data, and verified the results via the model- 
independent routine "Devin" . The approach is described 
below in more detail. 

To obtain the complex dielectric function, we fitted 
our transient reflectivity data to a highly flexible Drude- 
Lorentz function [S5]; the underlying Drude-Lorentz 
model of the optical data is based on the knowledge of 
the static optical spectra between a few meV and 4.6 eV, 
data taken from Ref . [S6] . 

For the fitting of the transient data, the Lorentz oscil- 
lators in our experimental range are allowed to change in 
order to reproduce the dynamical reflectivity. This pro- 
cedure yields the transient t\ and €2, displayed together 
with the static spectra in Fig. S2. In our experimental 
frequency range, the real and imaginary part of the di- 
electric constant have rather similar absolute values (Fig. 



10 



S2^4 and E). For this reason, the optical reflectivity, de- 
fined as R =| | 2 , is equally sensitive to the reactive 
and the absorptive components of the dielectric function. 
As a consequence, the complex e obtained through the 
Drude-Lorentz fitting of the optical reflectivity matches 
very well the one obtained via spectroscopic ellipsometry 
(see Fig. S2^4 and E and Ref. [S6] ) . The transient dielec- 
tric function is displayed in Fig. S2 (B-D, F-H) for light 
polarized along the a-axis and the c-axis at 10 K, and for 
light polarized along the c-axis at room temperature. 

Such a Drude-Lorentz modeling of the data also yields 
the spectral weight; the overall static N e ff spectrum for 
both light polarizations is displayed in Fig. S3. In the 
same plot, we show the derivative of N e f f as a function of 
the number of carriers which exhibits an inflection point 
between 1 and 1.5 eV for light polarized along the c-axis; 
this inflection point is found to be between 10 and 15 
carriers. Since Lu has a bd l electronic configuration, and 
we found a strong hybridization with at least one or two 
Ir 5d bands, we expect between 10 and 20 electrons to 
form the material conduction band and therefore the low 
energy c-axis conductivity should be mainly sensitive to 
these orbitals. 

By integrating the transient conductivity, we calcu- 
lated N e ff at each time step and obtained the time evo- 
lution of N e f f via Drude-Lorentz model, presented as 3D 
color maps in Fig. 2(A-C) of the main article. 

In order to verify that the transient changes of the SW 
we obtained by extrapolating the Drude-Lorentz model 
down to low frequency are robust and do not depend on 
the latter extrapolation, we fed into the "Devin" routine 
the real and imaginary parts of the dielectric function 
as an input in order to generate the SW up to a given 
frequency cut-off. Through this procedure it is also pos- 
sible to obtain an estimate of the uncertainty related to 
the knowledge of the spectrum being in a limited range 
with a limited resolution and with a given noise level. Be- 
cause we are interested in the relative changes of the SW 
induced by pump pulses, we use as an error our estimated 
accuracy in the determination of the transient dielectric 
function, which is of the order of 10 -4 , while the reso- 
lution of our spectrometer is approximately 1 meV. The 
transient dielectric function is fed into "Devin" and the 
model-independent SW changes are shown in Fig. 2(A- 
B) of the main article as white lines with errorbars. The 
fact that both these analyses give consistent results vali- 
date our description of the transient optical reflectivity. 

The reason why we can successfully obtain informa- 
tion on the CDW-induced spectral weight transfer from 
a limited spectrum in the visible light region is that, in 
this sample, the optical anisotropy, related to the Lu 
one-dimensional chains, strongly manifests itself around 
2 eV. This is visible in the static reflectivity spectrum 
(Fig. 1 of the main article), and also in the temperature- 
dependent dielectric constants obtained via ellipsometry 



and reported in Ref. [S6], which in fact show a partic- 
ularly large CDW-induced change between 1.5 and 2.5 
eV. 

Also, in strongly correlated solids, a phase transition 
can induce spectral weight transfers across a very broad 
energy range, extending to the visible light region. In 
these cases, it has been shown that the spectral region 
in the proximity of the screened plasma frequency of the 
material is very sensitive to these shifts [S7, S8]. Here, 
the screened plasma frequency along the c-axis (direction 
of the charge-density-wave) lies within our experimental 
range (zero crossing of ei in Fig. S3j4). 

THREE-TEMPERATURE MODEL 

The two-temperature model was first introduced by 
Kaganov et al. [S9] in order to describe the energy trans- 
fer between the electron and lattice subsystems, in the 
case where the electrons are driven to a higher tem- 
perature than the phonons. Its validity relies on the 
electronic temperature being of the order of the De- 
bye temperature (366 K in Lu5lr 4 Siio [S10]) and on the 
isotropy of the electron-phonon (e-ph) coupling func- 
tion. For anisotropic materials like cuprates [Sll] or 
iron-pnictides [SI 2], a selective coupling between elec- 
trons and a subset of the total phonon modes may be 
taken into account using the three-temperature model 
(3TM), governed by the following equations: 

2C e ^ = «/(*) + - 9(T e T h ) 

aC L fi = g(T e - T h ) - g c (T h - T c ) (3) 
(l-a)C L ^f = g c (T h -T c ) 

where T e , Th and T c are the temperatures of electrons, 
efficiently coupled phonons ( "hot phonons" ) and remain- 
ing modes ("cold phonons"), respectively. C e =jT e is the 
electronic specific heat, Cx is that of the lattice, a is 
the fraction of efficiently coupled modes, R is the static 
reflectivity and l s is the penetration depth (both at the 
pump energy). I(t) the pumping intensity, and K e is the 
electronic thermal conductivity (which is ignored in our 
case). Indeed, the electronic diffusion term ^(^e^f) 
changes neither the relaxation rates nor the e-ph thermal- 
ization time. The constant g governs the energy transfer 
rate from electrons to hot phonons, and is related to the 
second moment of the Eliashberg function A (w 2 ) through 
g = ^§JA (a; 2 ) [S13], A being the dimensionless e-ph cou- 
pling constant. g c governs the energy relaxation from hot 
phonons to the remaining lattice. 

To perform the 3TM simulations, we used an itera- 
tion procedure, calculating at each time step the depth- 
dependent temperature profiles. At each depth and time 
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step, we iterate the electronic and lattice part of the spe- 
cific heat. The electronic part is taken as C e = jT e , with 
7=23.42 mJ ■ mol^ 1 ■ K~ 2 [S10], and the lattice part is 
taken from Ref. [S14]. 

The e-ph coupling constant A was obtained via these 
simulations, given a logarithmic average of the phonon 
energies (w 2 ) of 30 meV, obtained by static optical spec- 
troscopy in Ref. [S6] . The e-ph coupling constant should 
in any case depend on the probing energy; as a check, we 
simulated the experimental data for several probe wave- 
lengths in order to show the result constancy in the ex- 
perimental energy range. The constant obtained through 
these simulations is given in Fig. S4C, from which we 
may notice its energy-independence in the range 1.6 -2 
eV. 

In fact, as this energy range is far from the electronic 
inter-band transitions (except very close to 2 eV, see the 
Lorentz oscillator positions in Fig. IB of the article), the 
probe is mainly sensitive to intra-band transitions. 

Considering N e ff, both experimental and theoretical 
determinations show that the carriers involved in the op- 
tical conductivity at 1.5 eV are mainly the Lu 5d and 
some Ir 5d, which mainly belong to the conduction band 
of the material; in fact, the plasma edge of the solid is 
close to 1.7 eV along the c-axis. Therefore, even though 
the numerous possible inter-band transitions should give 
a non-negligible contribution to the optical absorption, 
screening effects prevent the existence of long lived ex- 
cited states at those energies (1.5 eV) in such a metallic 
system, and the typical electron-hole lifetime should not 
exceed 10-20 fs. Therefore, on the time-scales of our ex- 
periments, > 50 fs, most if not all of the light excitation 
creates electron-hole pairs which relax instantaneously 
back in the conduction band with an excess energy that 
can be described by an effective temperature. Conse- 
quently, applying the 3TM in this energy range allows 
the measurement of the coupling between electrons at 
the Fermi level and lattice. 

The linear dependence of the fast relaxation time as a 
function of the maximum electronic temperature reached 
(see Fig. S4A), as expected considering this relaxation as 
due to e-ph coupling [S13], reflects also the applicability 
of the 3TM in our case. 

We may notice the excellent agreement between the pa- 
rameters obtained from this effective temperature model 
and the calculated partial e-ph coupling constants pre- 
sented in the following section. Indeed, both the fraction 
of efficiently coupled phonons (Fig. S4D) and the sum of 
their partial e-ph coupling (Fig. S4C) may be retrieved 
by 3TM simulations, enhancing our confidence in the re- 
liability of this model in describing our transient data. 



ELECTRONIC STRUCTURE CALCULATIONS 

Two separate sets of calculations were performed, the 
first using the highly accurate all-electron full-potential 
linearized augmented plane-wave method (often referred 
to as FP-LAPW), as implemented in the ELK code [S15], 
and the second using pscudopotcntials, as implemented 
within the QUANTUM ESPRESSO code [S16], for the 
calculation of the phonons and e-ph coupling. 

The lattice constants used were the experimental ones 
(a = 12. 4936 A, c = 4. 1852 A), as were the atomic po- 
sitions [S17], giving a unit cell of two formula units (38 
atoms). The calculations made using the ELK code [S15] 
had a cutoff (for planewavcs in the interstitial region) de- 
termined by fc max = 7.0/i? m in, where i? m i n is the smallest 
muffin-tin radius. These muffin-tin radii were Lu (2.3328 
a.u.), Ir (2.3654 a.u.) and Si (1.8369 a.u.). The self- 
consistent cycle was carried out on a mesh of 64 fc-points 
in the irreducible (one sixteenth) of the Brillouin zone, 
and the Fermi surface evaluated from the converged po- 
tential on a denser mesh of 1936 fc-points within one 
eighth of the Brillouin zone. 

The calculations presented are scalar-relativistic and 
omit spin-orbit coupling (after the effect of its inclusion 
was investigated and found to be negligibly small at the 
Fermi level.) The exchange-correlation functional was 
that of Perdew- Wang/ Ceperley- Alder, Ref. [SI 8]. 

The calculated band structure in the energy range [-8 
eV, +2 eV] is given in Fig. S5A. In this range, the latter 
is composed of strongly hybridized Lu 5d and 4/, Ir 5d 
and Si 2p, the Lu 4/ giving localized density of states 
around -4.5 eV, while electrons at the Fermi surface come 
mainly from Lu 5d, Ir 5d and Si 2p as shown in Fig. 
S5B and Fig. 5 of the main article. A large amount 
of "interstitial" electronic density is also found at the 
Fermi level, indicating their delocalized nature as well as 
the strong hybridization degree in this system. 

The e-ph coupling resulting from the mass enhance- 
ment may be calculated from the density of states (DOS) 
at the Fermi level, following the formula N(Ep) = 
N b (E F )^ = N b (E F )(l + A), where N(E F ) is the renor- 
malizcd DOS and Nb(Ep) the bare one, both at the Fermi 
level E F ', m is the free electron mass and m* the effective 
one. The renormalizcd DOS obtained by specific heat 
measurement being 0.26/eV/atom/spin, (i.e. 19.76/eV), 
and the bare DOS calculated here of 11.02/eV if one con- 
siders that 30% of the Fermi surface is gapped, one ob- 
tains (1 + A) = 1.79 giving an e-ph coupling constant 
A = 0.79. 

We also present in Fig. S6 the real and imaginary 
parts of the calculated susceptibility for two bands, cor- 
responding of the Fermi surfaces indicated by arrows in 
Fig. 5^4 of the main article. One may notice the peak 
at two different wavevector, in both real and imaginary 
parts. These peaks appear at nesting wavevectors which 



12 



are different between the two presented bands, indicating 
the formation of several charge-density waves. 

For the phonon and e-ph coupling calculations, lin- 
ear response within the pseudopotential approach pro- 
vided by the QUANTUM ESPRESSO code was used 
[S16]. The Lu pseudopotential had the / electrons as 
core states, and the resulting electronic structure was 
verified against that produced by the all-electron ELK 
calculation. Both ultrasoft and norm-conserving pseu- 
dopotentials were used, the latter needed for the calcula- 
tion of the dielectric function with the EPSILON package 
within QUANTUM ESPRESSO. The cut-offs for the en- 
ergy and charge density were 40 Ry and 400 Ry, respec- 
tively. Phonons at the T point were calculated from self- 
consistent calculations performed on a (6,6,16) fc-point 
mesh and e-ph parameters were evaluated on a denser 
mesh of (12,12,32). For the calculations of dielectric 
function, the routine EPSILON within the QUANTUM 
ESPRESSO package was used. Here, Brillouin zone inte- 
grations were carried out on a mesh of over 6000 fc-points 
with a Gaussian smearing of 0.25 eV applied to the eval- 
uation of the inter-band contribution. 

We present the partial e-ph coupling constants for each 
of the 114 modes of the system as a function of the mode 
energy at q = in Fig. S7B. Two strongly coupled modes 
are found having respectively a partial e-ph coupling con- 
stant of 0.166 and 0.131, involving mainly Si atom motion 
since the latter are the lighest elements of the crystal 
structure. One of these two modes is represented with 
yellow arrows in Fig. 5C of the main paper. The sum 
over all modes gives A = 0.97, so in good agreement with 
the value obtained from the DOS renormalization at the 
Fermi level (when one considers that the electron-phonon 
calculation has only been made at the T point). 

The histogram of partial e-ph coupling constants is 
shown in Fig. S7A. From such a representation, one is 
able to extract the partial e-ph coupling of some per- 
centage of the most efficiently coupled modes, as repre- 
sented with a red-dashed line in Fig. S7A. This area con- 
tains 21 modes, (or 18 % of the phonons), and summing 
up their A constant one obtain 0.62. It indicates that 
the 18 % of the phonon modes having the largest e-ph 
coupling have, all together, an e-ph coupling constant of 
0.62. This observation is in excellent agreement with the 
results of 3TM simulations, also indicating the reliabil- 
ity of these electronic structure and phonon calculations. 
This value of A gives a superconducting transition tem- 
perature of ~ 9 K using the McMillan formula [S19], as 
experimentally observed in these system above the crit- 
ical pressure, and suggests that conventional e-ph medi- 
ated superconductivity takes place in Lush^Siio, obey- 
ing Bardeen-Cooper-Shrieffer theory. Indeed, since upon 
photoexcitation the CDW is molten, and therefore does 
not compete anymore with superconductivity (as in the 
system under pressure), obtaining a similar A value may 
reflect effectively the order parameter giving rise to su- 



perconductivity and not the one related to Peierls order- 
ing. 
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FIG. 7: (A) and (E) Static dielectric function for both light polarizations. Color maps: transient dielectric function at different 
temperatures and along the different orientations. 
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FIG. 9: Three temperature model analysis. (A) Fluence dependence of the signal obtained at 10K, with pump and probe 
parallel to c-axis, at an energy of 1.7 eV; in the inset the fast decay time, associated with the electron-phonon coupling time, 
is plotted as a function of the maximum electronic temperature reached in average in the skin-depth together with its linear 
fit. (B) Electron, hot phonon and cold phonon temperatures as a function of time delay obtained by simulating the transient 
reflectivity shown in open circles in the inset, with the simulation as the solid line; this curve corresponds to a geometry where 
both pump and probe were parallel to the c-axis, at a pump fluence of 3.1 mj/cm 2 , and a probing energy of 1.7 eV. (C) 
electron-phonon coupling constant obtained from the simulations, for probe polarization parallel to c-axis or within the (a, b) 
plane, simulated at a pump fluence of 3.1 mj/cm 2 as a function of the probing energy; (D) corresponding fraction of coupled 
modes. 



15 




FIG. 10: Electronic band structure. (A) Calculated electronic band structure using the ELK code. (B) Density of states (note 
that the contribution of the localised Lu 4/ states is truncated). 




q (in 2rc/c units) 

FIG. 11: Calculated real and imaginary part of the susceptibility, for two of the bands close to the Fermi level (labeled 278 and 
280). 
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FIG. 12: Calculated electron-phonon coupling. (A) Histogram of the electron-phonon coupling constants A; the red area 
represents the 18% more coupled modes, whose total lambda gives 0.62. (B) Energy distribution of the A. 



